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Ground state of the one-dimensional transverse field Ising model is investigated under the 
hyperbolic deformation, where the energy scale of j-ih bond is proportional to the function 
cosh [j A] that contains a parameter A. Although the Hamiltonian is position dependent, the 
ground state is nearly uniform and finitely correlated. We observe the energy cross over between 
the ordered and disordered state with respect to the transverse field. The model shows first 
order phase transition, and the discontinuities in the magnetization and entanglement entropy 
at the transition point detect the Ising universality. 
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1. Introduction 

Spacial uniformity is a fundamental concept in physics. 
If a Hamiltonian of a system is uniform, in the manner 
that it is represented as a spacial sum or integral of posi- 
tion independent local terms, the corresponding ground 
state is expected to be uniform, provided that there is 
no spontaneous symmetry breaking that causes spacial 
modulation. How about the opposite? When the ground 
state is uniform, is it expected that the Hamiltonian is 
also uniform? This is not true; a nonuniform Hamiltonian 
can possess a uniform ground state. A trivial example is 
the system that consists of independent spins. Consider 
the position dependent spin Hamiltonian 



(i.i) 



where T > represents the external magnetic field to 
the ^-direction, a site dependent positive factor at j- 
th site, and the Pauli operator. The corresponding 
ground state is the complete ferromagnetic state, where 
all the spins are polarized to the x-direction. 

The example in Eq. (1.1) might be too trivial, since 
there is no inter-site couplings. Thus let us consider the 
transverse field Ising (TFI) model on the one-dimensional 
(ID) lattice, as a more realistic reference system. We 
treat the position dependent TFI model defined by the 
Hamiltonian 







(1.2) 



where J > represents the longitudinal nearest-neighbor 
coupling, and / • is a site dependent positive factor. When 
all the factors fj and g • are equal to unity, the ground 
state of this system is uniform in the thermodynamic 
limit. In this uniform case there is a quantum phase tran- 
sition at r = J, where the model becomes self dual. 1 ) 
(See appendix.) 

In this article we focus on the case where the posi- 
tion dependence is given by / • = cosh [j A] and g- = 



cosh [(j — I) A] , where A is a nonnegative parameter. 
The corresponding Hamiltonian 



j 

-r£coBh[(j-|)A] a 



(1.3) 



for the case A > can be interpreted as the one parame- 
ter deformation — the hyperbolic deformation — to the 
uniform case when A = 0. The 'bulk' part of the ground 
state of H C (X) is expected to be uniform even when 
A > 0, since the path-integral representation of the imag- 
inary time evolution by H C (X) would be given by uniform 
classical action on the 1 + 1 hyperbolic plane. 2 ) Such a 
uniformity under hyperbolic deformation has been ob- 
served for the deformed 5 = 1/2 and 5 = 1 Hciscnbcrg 
spin chains. 3 ' 4 - 1 We confirm this uniformity for the case 
of the deformed TFI model in the next section. 

We employ the density matrix renormalization group 
(DMRG) method, 5 ~ 8 ' and obtain the lowest energy states 
of finite size systems under both ferromagnetic and para- 
magnetic boundary conditions. A natural interest in the 
ground state of H C (X) is the ordered-disordered transi- 
tion. A classical analogue of the deformed TFI model is 
the classical Ising model on the two-dimensional (2D) 
hyperbolic lattices, 9 ) which exhibit the mean-field like 
second-order phase transition. 10 ~ 21 ) How about the de- 
formed TFI model? We investigate the spontaneous mag- 
netization in §3. In contrast to the classical cases, the 
deformed TFI model exhibits first-order transition. The 
A-dependence of the quantum entropy is also observed. 
Conclusions are summarized in the last section. 

2. Uniformity in the Ground State 

Under the hyperbolic deformation, the energy scale in 
the Hamiltonian H C (X) blows up exponentially with \j\. 
In order to well define the eigenvalue problem for H C (X), 
we consider finite-size systems. For simplicity, we treat 
the cases where system size L is even, and label sites from 
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j = —L/2 + 1 to j = L/2. Therefore the center of the 
system is between j = and j = 1, where the strength 
of the nearest neighbor interaction is the smallest. 

Boundary condition is essential for the determination 
of the ground state, since the ratio of the boundary (or 
surface) energy with respect to the total energy does not 
vanish in the large L limit. 11 ' 12 ) We have to aware of 
this characteristic behavior of the hyperbolic deforma- 
tion when we consider the ground-state phase transition. 
We choose either the paramagnetic boundary condition 
or the ferromagnetic one. The former is imposed by fixing 
the spins of both ends to the x-direction. Since the inter- 
site coupling in the TFI model is mediated only by lon- 
gitudinal interaction <7J<tJ +1 , the paramagnetic bound- 
ary condition decouples the boundary sites at the both 
ends j = —L/2 + 1 and j = L/2 from the inner part 
—L/2 + 2 < j < L/2 — 1. As a result, the effect of param- 
agnetic boundary condition is represented by the Hamil- 
tonian 



#P;l(A) 



j=L/2-2 

-J ]T cosh[jA] a]a) +x (2.1) 

j=-L/2+2 



j=L/2-l 

£ 

] = -L/2+2 



COsh[(j - \) A] (Tj 



where there is no constraint for all the spins that appear 
in the above equation. On the other hand, the ferromag- 
netic boundary condition is imposed by fixing the bound- 
ary sites to the ^-direction. In this case the Hamiltonian 
is expressed as 



#f : l(A) 



#p ; l(A) 



(2.2) 



+ COsh [(f - !) A ] (V-L/2+2 + °l/2-l) > 

which is equivalent to put magnetic field to z-direction 
only at the position j = —L/2 + 2 and j = L/2 — 1. 
Thus under both boundary conditions we effectively treat 
(L — 2)-site system in numerical analysis. We choose the 
parameter J as the unit of energy throughout this article. 

We employ DMRG method 5-8 - 1 for the numerical deter- 
mination of the ground state. A direct application of the 
finite-system DMRG algorithm encounters a numerical 
instability, which is caused by the blow-up of the energy 
scale cosh [j A] with respect to \j\. In order to stabilize 
the computation, we treat 

H§,. L {\) = H$. L (\) - (m. L (X)) 

^f ; l(A) = H$. L {\) - (H$. L (\)) (2.3) 

instead of Hp. L (X) and Hp. L (X) directly, where ( ) de- 
notes expectation value taken by the ground state. The 
smallest eigenvalue of both Hp. L (X) and Hp. L (X) is zero 
by definition. Since the ground state is not a priori 
known, the subtraction process in Eq. (2.3) is performed 
iteratively. The decay of the density matrix eigenvalue is 
rapid, 21 ) and for all the range of deformation parameter 
we examined 1/32 < A < 1 it is sufficient to keep to = 16 
block-spin states; even for the worst case A = 1/32 the 
truncation error is of the order of 10~ 10 in the density 
matrix eigenvalues. 



V 



X=1,L=32, m=16 
r=0.96 



r=i.04 



Fig. 1. Expectation value of h- - +l in Eq. (2.4) for the ground 
state when Y = 0.96 and Y = L04. 



For the analysis of of local energy, we employ a bond 
operator 



h 



2 a > 



J(7 i cr j+i 



r 

r 



(2.4) 



which coincides with the local energy of the j-th bond 
of the undeformed TFI model. Figure 1 shows the ex- 
pectation value (hj j +1 ) calculated for the ground state 
when A = 1 and L = 32. We choose two typical cases, 
r = 0.96 with ferromagnetic boundary condition and 
T = 1.04 with paramagnetic one. Despite of the very 
strong position dependence in the Hamiltonian, bound- 
ary correction decays rapidly, and (hj j +1 ) is nearly uni- 
form deep inside the system. Figure 2 shows the on-site 
transverse interaction (IVJ), which is nearly uniform in- 
side. We also observed similar uniformity for other local 
observables such as (crj) and (ctJ<tJ +1 ). The rapid decay 
of the boundary corrections is observed for all the values 
of the transverse field T we have examined. The behavior 
suggests that the ground state of the hyperbolically de- 
formed TFI model is always finitely correlated, and that 
the bulk part of the ground state is well approximated 
by the uniform matrix product. 22, 23 ) 



X=1,L=32, m=16 
r=0.96 



-15 -10 -5 5 10 15 
j 

Fig. 2. On-site transverse interaction (Ycr?). 



3. Phase Transition 

In order to capture the nature of ground-state phase 
transition, let us compare the energy of the ferromag- 
netic state with that of paramagnetic one. It is, however, 
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difficult to directly compare (Hp. L (X)) with (_ffp\ L (A)), 
since subtraction of the boundary energy is not straight- 
forward. Therefore we use the expectation value (/io,i) 
at the center of the system as a representative value for 
the local energy density of the bulk. The uniformity ob- 
served in the previous section would justify this way of 
evaluation of the ground-state energy. Figure 3 shows 
(h a _i) for L — 8, 16, 24, and 32 when A = 1. Solid lines 
and dotted lines, respectively, represent data calculated 
with ferromagnetic and paramagnetic boundary condi- 
tions. System size dependence between L = 24 and 32 is 
negligible, and for these cases the solid and dotted lines 
crosse at T = 1. Although the deformed Hamiltonian is 
not invariant under the duality transformation, still the 
lowest energy state alternates at the point r = 1. (See 
Appendix.) 





-1.23 






-1 .24 1 






-1.25 




A 








-1.26 




o" 








-1.27 




V 








-1.28 


X=1,m=16 '*v>./^^ 






L=8 + 




-1.29 


/. 15 * ''--.'."} 




L=24 □ 






L=32 




-1.3 





0.96 0.97 0.98 0.99 1 1.01 1.02 1.03 1.04 



Fig. 3. Energy density {h x ) at the center of the system. Solid 
lines and dotted lines, respectively, represent data calculated 
with ferromagnetic and paramagnetic boundary conditions. The 
L-dependence is very small for L > 16. 



It should be noted that in a certain region of T around 
T = 1 the ferromagnetic and the paramagnetic states co- 
exist even in the large L limit, and one of these states 
is chosen by the imposed boundary condition. Thus one 
should be careful about the definition of the phase tran- 
sition. Since we have considered that (h ,i) represents 
the bulk property of the system, it is natural to choose 
the state that gives lowest value of (/io,i) as the 'stable' 
ground state. The ground-state energy thus determined 
has a kink at the transition point T = 1 . We checked the 
presence of kink in (/io,i) for the cases of weaker deforma- 
tion down to A = 1/32, where the sufficient system size L 
to obtain the convergence in (fto,i) grows roughly propor- 
tional to 1 /A. The fact suggests that the hyperbolic defor- 
mation introduces a length scale £ oc 1/ A to the system, 
which coincides the radius of curvature R oc — 1/A of the 
corresponding hyperbolic plane. 2,9 ) Actually, the length 
£ already appeared as the dumping length in Fig. 2. 

Observation of the order parameter clarifies the na- 
ture of phase transition. Figure 4 shows the magneti- 
zation (ctq) at the center of the system j — 0, plotted 
for the states that gives lowest value of (h 01 ); the fer- 
romagnetic boundary condition is imposed when T < 1, 
and the paramagnetic condition when T > 1. We show 
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Fig. 4. Order parameter (eg) observed at the center of the sys- 
tem. 



the results when A = 1, 1/2, 1/4, 1/8, and 1/16 calcu- 
lated for sufficiently large system size L. There is a finite 
jump in (<7q) at the transition point T = 1, and thus the 
phase transition is first order. The A dependence of this 
jump is plotted in Fig. 5. The eighth power of the jump 
is proportional to A. This dependence is consistent with 
the Ising universality. The length scale £ cx 1/A intro- 
duces an effective deviation to the transverse field of the 
amount 



Ar cx = j oc A 



(3.1) 



from the criticality of the uniform TFI model, where 
v = 1 is the critical exponent for the correlation length. 
The observed jump in the spontaneous magnetization is 
therefore proportional to (Ar) 1 / 8 oc A 1 / 8 . 
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Fig. 5. Jump in the magnetization at T = 1 with respect to A. 



The bipartite entanglement entropy S provides sup- 
plemental information for the phase transition. Figure 6 
shows the entropy measured at the center of the system. 
Consistent with the spontaneous magnetization, the en- 
tropy does not show singular behavior around the tran- 
sition point. Figure 7 shows the entropy at T = 1 for 
the ordered state (A) and that for the disordered state 
(B), and also their sum (A)+(B) and difference (B)-(A). 
We draw a fitting line —(1/6) log A + 0.2475 for the sum 
(A)+(B). This suggests that the average of the entangle- 



4 J. Phys. Soc. Jpn. 



Full Paper 



H. Ueda, A. Gendiar, V. Zauner, T. Iharagi, and T. Nishini 



0.6 



o 

H 

a 
W 

a 



g 

S3 

I 



0.2 



0.1 



X=1/2 
- X=1/4 
X=1/8 
X=1/16 


1 


° • : 


o • 


- m=16 






x x x x 


i 8 1 = 


+ + + + 


+ + + 




.9 


0.95 


1 

r 


1.05 1 



Fig. 6. Bipartite entanglement entropy measured at the center of 
the system. 



ment entropy [(A)- 



(B)]/2 at T = 1 is expressed as 
0.2475 



1 1 1 

2'6 l0g A + 



(3.2) 



for a wide range of A. This dependence coincides with the 
fact that the leading term of the entanglement entropy 
is expressed as (c/6) log W, where W is the system size 
and c is the central charge. 24 ) This is because 1/A is 
proportional to the length scale £, and the TFI model 
belongs to the class where c = 1/2. 
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Fig. 7. Entanglement entropy with respect to A. 



The fitting to the difference (B)-(A) gives the estimate 
0.3426 in the small A limit. This value is about the half 
of log 2 ~ 0.6931. The result suggests that difference of 
entanglement entropy between ordered and disordered 
states is a constant, if the correlation lengths of the both 
states are the same. We conjecture that the difference 
~ (1/2) log 2 captures the boundary condition, but the 
detail is not clarified yet. 

4. Conclusion and Discussion 

We have introduced the hyperbolic deformation to the 
ID TFI model. It is shown that inner part of the ground 
state is uniform, as was observed other systems under 
hyperbolic deformation. 2 ~ 4 ) The deformation introduces 
a characteristic length £, which is proportional to 1/A, 
and the ground-state phase transition becomes 1st order. 
Discontinuities in both spontaneous magnetization and 
bipartite entanglement entropy are consistent with the 



Ising universality class. 

The observed 1st order transition is different from the 
mean-field like 2nd order transition observed for the clas- 
sical Ising model on hyperbolic lattices. 10 ~ 21 ) Thus spa- 
tial anisotropy may change the nature of the phase tran- 
sition. To clarify the difference, we have to construct a 
well defined quantum-classical correspondence by some 
means such as the Trotter decomposition. 25,26 ' 
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Appendix: Duality Relation 

The undeformed TFI model has a symmetry mediated 
by the duality relation. ^ Let us introduce a new set of 
Pauli operators Tj , rj , and rj . Then the Pauli operators 
crj, ctJ, and <rj can be expressed as 

a 3 = T /-1 T / 




(AT) 



where the transformation is nonlocal. One can verify the 
relation <7jcrJ +1 = rj. Substituting the above relation to 
H C (A) in Eq. (1.3), we obtain the Hamiltonian 



IF (A) = -J ^ cosh [j A] 

3 

-r£cosh[(,-i)A] t^t/ (A.2) 

3 

written by the new set of Pauli operators. When A = 
the transformed Hamiltonian has the same form as the 
original Hamiltonian, where parameters J and T are ex- 
changed. For finite A, the transformed Hamiltonian in 
Eq. (A-2) has the form where lattice indices in the defor- 
mation function is shifted by 1/2. The observed energy 
crossover at T — 1 in Fig. 3 suggest that this shift is not 
essential for the bulk part of the ground state. If we shift 
the index by 1/4 in advance and define the Hamiltonian 
as 



IF' (A) = -J^cosh[(i + i)A] a]a] +x 

3 

-r£[cosh(j-i)A] o*, 

3 

then the transformed Hamiltonian is given by 
IF' (A) = -J^cosh[(i + i)A] r? 



(A-3) 
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-r^cosh[(i-I)A] 

3 

One finds the self duality between Eq. (A-3) and 
Eq. (A-4) when J = T even when A is positive. 
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